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Abstract 


The  problem  considered  here  is  that  of  the  transition  region 
separating  a  uniform  Maxwellian  Plasma  from  its  confining  magnetic 
field.   The  geometry  chosen  is  one-dimensional.   The  magnetic  field 
profile  within  the  boundary  layer  is  determined  by  the  current  dis- 
tribution, i.e.,  by  the  paths  of  the  ions  and  electrons.   The  paths 
of  these  particles  on  the  other  hand  are  determined  by  the  fields  in 
which  they  move.   The  equation  of  motion  of  the  ions  is  treated 
exactly  and  the  electrons  in  the  guiding  center  approximation.  A 
self  consistent  solution  of  Maxwell's  equations  is  obtained.   These 
results  are  compared  with  those  obtained  for  a  mono-energetic  dis- 
tribution and  found  to  agree  to  within  10^. 


1.  Introduction 

This  report  is  an  extension  of  an  earlier  work  of  one  of  the  authors'. 
The  problem  discussed  there  is  that  of  the  transition  region  separating 
a  uniform  plasma  and  its  confining  magnetic  field.   It  was  possible  to 
perform  the  calculations  exactly  due  to  the  simplifying  restriction  that 
the  plasma  be  mono-energetic,  i.e.,  that  all  particles  have  the  same  energy. 
If  both  the  electron  and  ion  component  are  treated  exactly  the  analysis 
is  quite  laborious.   On  the  other  hand  if  the  ions  are  treated  exactly  and 

the  electrons  in  the  guiding  center  approximation  the  error  involved  in 

6 
the  transition  layer  thickness  is  only  5  parts  in  10  . 

In  the  present  paper  we  should  like  to  consider  a  more  realistic 
plasma--a  plasma  which  becomes  Maxwellian  far  from  the  transition  region. 
Unfortunately  it  is  not  possible  to  solve  this  problem  exactly.   However 
the  agreement  between  the  exact  and  approximate  calculation  for  the  mono- 
energetic  plasma  suggests  that  the  approximate  calculations  should  suffice 
for  the  Maxwellian  plasma  as  well. 

We  shall  not  dwell  here  on  the  formulation  of  the  problem  or  upon  the 
nature  of  the  approximation  as  these  points  are  discussed  in  detail  in  I. 
We  shall  content  ourselves  with  a  brief  statement  of  the  problem  and  a 
presentation  of  the  results. 

2 .  Statement  of  the  Problem 

All  quantities  are  assumed  to  depend  only  on  a  single  variable  x. 


J.  Hurley,  Physics  of  Fluids  6,  83  (1963) .   Hereafter  referred  to  as  I. 


The  magnetic  field  is  directed  everywhere  in  the  z  direction  and  is  uniform 
from  X  =-  00  to  X  =  0.  As  we  cross  the  plane  x  =  0  we  enter  into  the  domain 
of  the  collision  free  plasma  where  the  magnetic  field  decreases  and  eventual- 
ly becomes  uniform  again  at  x  =  +  oo  where  the  plasma  is  required  to  be 
Maxwellian.  We  further  require  that  the  ion  density  be  identically  equal 
to  the  electron  density  throughout  the  plasma  domain.   (Thus  there  are  no 
electric  fields.)   The  paths  of  the  particles  are  determined  by  the  magnetic 
field  in  which  they  move.  The  field  in  turn  is  determined  by  the  paths  of 
the  particles.   The  problem  is  to  solve,  self-consistently,  the  equations 
of  motion  of  the  particles  (Liouville's  equation)  and  Maxwell's  equations. 

3.   Analysis 

The  magnetic  field  is  determined  by  the  relation 

dB     -Un   .       ,  ,  . 

—  =  (.j+  +  J.)  -'        (1) 

dx      c 

where  j   and  j   are  the  ion  and  electron  currents  and  are  given  by 

j^  =  -e  /  y  f^  (x,  X,  y,z)    di  dy  dz  (2) 

and  f  and  f  are  the  ion  and  electron  distribution  functions.  The 
distribution  functions  will  satisfy  Liouville's  equation  if  they  are 
chosen  to  be  functions  of  the  constants  of  the  motion. 


These  constants  are 


P+   =  m^y  +  e  A/c  (3) 


where  A(x)  is  the  vector  potential  defined  by  the  relation  B  =  dA/dx 
and  normalized  so  that  A(o)  =  0.  We  write  f(E,P,Q)  for  f  (x,x,y,z) 
when  expressed  as  a  function  of  E,  P,  and  Q.  We  shall  choose  as  the 
ion  distribution  function 

f^   =   (rny23tkT)  "  N^(oc)  T]    (P^  -  2m^E)  e"  ^/^^      (h) 

where  T]  is  the  unit  step  function  and  N  (<»)  is  the  ion  density  at 
X  =  +  00  (as  may  be  seen  by  integrating  f  ) .   The  variables  in  (k)    are 
restricted  by  the  inequality 


2  m  E  -    {F   -  —   f   -   q^  >  0 


It  may  be  seen  from  (2)  that  this  is  equivalent  to 


Such  a  choice  for  the  ion  distribution  function  requires  (l)  that  there 
shall  be  no  particles  in  the  domain  x<0  (or  equivalently  in  the  domain 


A  <  0  since  we  assume  that  A(x)  is  a  monotonic  function  of  x)  and 
(2)  that  the  distribution  shall  approach  a  Maxwellian  distribution  as 
X  approaches  infinity. 

We  may  now  compute  the  ion  current  required  to  solve  the  differ- 
ential equation  (l) .   Substituting  (k)    into  (2)  we  find 

j^   =    (e/m)    N_^(cx>)  J 2  m_^kT  a    (    1   -   Prob  a)  (5) 

where 

a  =  eA/c    /Bm^kT  (6) 


^        ^,  2  r  -U^        ^ 

Prob  a     =  /      e  du 


However  to  compute  the  electron  current  it  is  necessary  to  find  f_, 
the  electron  distribution  function.   This  may  be  done  by  solving  the 
integral  equation  obtained  by  equating  the  ion  and  electron  densities 


/ 


f_^  (x,  i,  y,  z)    di  dy  dz  (?: 


Unfortunately  this  is  too  difficult  to  solve  when  f_^  is  given  by  (1+) 
It  was  shown  in  I  that  it  is  possible  to  bypass  this  step  when  the 
guiding  center  approximation  is  used  for  the  electrons,  for  then 


-  5 


dvi'J  d  N  d  N, 


d  A 


k  T (8) 


where  p^"'^   is  the  Indicated  component  of  the  electron  stress  tensor. 
We  use  (h)    and  (7)  to  determine  N  and  substitute  into  (8)  giving 

ekT  W^  (00)      2 


yaTm_^kT 


(9) 


Introducing  (5)  and  (9)  into  (l)  we  obtain  the  ordinary  differential 
equation  for  the  magnetic  field 


^e 


L  a  (  1  -  Prob  a  )  +    e~^  /  2   J7~  ] 


(10) 


where 

i  =J2  n    e^  N^  (»)  /  m^c^   x  (ll) 

The  initial  conditions  A(o)  =  0  and  b(0)  =  B   become 


a  (0)  =  0   ,  (  ^  )^    =  B^  /  J   16  n  N^  (-)  kT         (l2) 
4=0 


We  shall  find  it  more  convenient  to  parametrize  the  plasma  by  its 
field  at  X  =  00  (b^)  rather  than  by  N_^  (00)  kT.  This  may  be  accom- 
plished by  integrating  (lO)  once  from  ^  =  0  to  4  =00.  After  some 
simplifications  we  obtain 


T.2    J2. 
B     B 

—2-   =  ^^  +  2  N   («.)  kT  (13) 

Srt    8rt 

which  we  recognize  as  the  condition  of  equilibrium  for  a  large  par- 
allelopiped  whose  left  face  lies  in  the  vacuum  at  x  =  0  where  the 
pressure  is  purely  magnetic  and  equal  to  B  /8jt   and  whose  right  face 
lies  in  the  Maxwellian  plasma  at  x  =  oo  where  the  field  pressure  is 
B  /  8jt  and  the  ion  pressure  W  (oo)kT  and  the  electron  pressure 


W  (oo)   kT  =  N_^(«>)  kT.   Ohserving  that  B  =  (dA/dx) 
write  the  second  initial  condition  (l2)  in  the  foi 


x=0 


we  may  re- 


r  a  Q:  n      _   / 


For  purposes  of  comparing  these  results  for  the  Maxwellian  plasma  with 
those  obtained  in  I  for  the  mono-energetic  plasma  we  have  integrated 
(lO)  -  (l2)  numerically  for  the  special  case  B  =  B  / /2~   .      The  re- 
sults are  shown  in  figures  (l)  and  (2).  As  a  measure  of  the  transition 
layer  thickness  we  chose  that  point  at  which  the  field  has  fallen  by 
•  9  of  the  total  discontinuity.  This  occurs  at  !|  =  1.15  or 


X  =  1.15  ,/m^c^/2ire^N_^  («>)   =  k.60   c  Jm_kT     /e   B^ 

where  we  have  made  use  of  (ll)  and  (13) .   Wow  the  average  energy  per 
particle  is  E  =  3/2  kT.   In  terms  of  this  energy  we  may  write  for  the 
transition  layer  thickness 


^o  =  3-77  c  Vm^E^  /  e  B^ 

This  is  to  be  compared  with  the  exact  solution  for  the  mono -energetic 
plasma 

^o  =  3.30  c  Jl^     I    e   B^ 

where  E_^   is  the  kinetic  energy  common  to  all  particles.  We  conclude 
then  that  the  depth  of  the  transition  layer  is  essentially  insensitive 
to  the  energy  distribution  in  the  plasma. 


\ 
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